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ABSTRACT 



The chaos in stellar systems is studied using the theory of dynamical systems and the Van Kampen stochastic differ- 
ential equation approach. The exponential instability (chaos) of spherical N-body gravitating systems, already known 
previously, is confirmed. The characteristic timescale of that instability is estimated confirming the collective relaxation 
time obtained by means of the Maupertuis principle. 
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1. Introduction 

Stellar dynamics, a well established discipline, continues to 
be an active field of research due to its complex dynamics. 
Newly developed sophisticated tools on nonlinear dynamics 
or ever increasing numerical capacities are readily applied 
to N-body gravitating systems. 

The problem of relaxation of gravitating systems is 
' among the key ones in stellar dynamics since it can be 
directly constrained by observations, especially, of well- 
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relaxed globular clusters and elliptical galaxies. 

Historically, plasma methods were among the first to 
be applied to gravitating systems (Chandrasekhar 1942), 
neglecting, however, a drastic difference between plasma 
and long-range gravity. In retrospect, it can seem strange 
how easily were ignored the challenges in theory and ob- 
servations, and later also in numerical simulations: (a) the 
Coulomb logarithmic cut-off and hence the canceling of the 
N-body effects was applied at the absence of Debye screen- 
ing; (b) the result contradicted observations, i.e. the plasma 
two-body relaxation timescale exceeds the age of elliptical 
galaxies by several orders of magnitude; (c) the two-body 
timescale could never be identified in numerical studies. 

Then comes the epoch of realization of importance of 
chaos, occasionally in unexpected forms, for nonlinear sys- 
tems. Chaos caused by small perturbations appears cru- 
cial even to nearly integrable problems, such as the dy- 
namics of planetary system, (see Laskar & Robutel 1993, 
Morbidelli 2002 and refs therein); for non- integrable N- 
body systems the situation is far more complex. The 
Fermi-Pasta-Ulam problem is another example of ap- 
parently simple, but still not well understood system 
(Fermi et al. 1955). 

Do chaos, chance and randomness have significant role 
also in the evolution of stellar systems, or the two-body 
plasma approach is the whole story? Although it has been 
generally agreed that chaos must affect N-body dynam- 
ics (see Gurzadyan & Pfenniger 1994, Contopoulos 2002, 
Regev 2008), proper treatments require well-founded ap- 



proaches 1 . Application of ergodic theory methods enabled 
us to prove a notable result: the spherical systems are expo- 
nentially unstable and that chaoticity drives their dynam- 
ics (Gurzadyan & Savvidy 1984, 1986). An exponential law 
defines an intrinsic timescale, the collective relaxation 
time, which for real stellar systems has a value interme- 
diate between the dynamical (crossing) and the two-body 
timescales. This was obtained by estimating the divergence 
of the trajectories of the system in a Riemannian space 
defined by the potential of the interaction (Maupertuis 
principle), i.e. by a method known in theory of dynam- 
ical systems (Arnold 1989). Importantly, the derived col- 
lective (N-body) relaxation timescale fits the observational 
data (Vesperini 1992). The three timescales, i.e., the dy- 
namical Tdyn, collective r cr , and two-body t& timescales, 
correspond to the 3 distance scales of the system, specif- 
ically, its size D, mean inter-particle distance d, and ra- 
dius of gravitational influence of particles r^, in particular 
(Gurzadyan & Savvidy 1984, 1986) (see also Lang 1999), 
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When the role of complex N-body dynamics was finally 
recognized, another confusion did appear, namely, in as- 
signing of the dynamical time as the relaxation timescale, 
i.e., for reaching fine-grained equilibrium. However, this 
again contradicts observations: globular clusters would then 
have already disappeared because of the evaporation of 
stars (Ambartsumian 1938) within 100 crossing times, i.e., 
within around 100 mln years. 2 The dynamical timescale 



1 Ruelle mentions the appearance of numerous incorrect pa- 
pers about chaos when chaos became a fashion (Ruelle 1993). 

2 One of the motivations for concluding that the dynami- 
cal time equals the relaxation time were the numerical experi- 
ments on the apparent divergence in real space first observed in 
(Miller 1964) for systems of between N=8 and 32 particles, for 
which obviously no relaxation process has any sense. Properly 
performed numerical experiments (El-Zant 1997) confirmed the 
timescale obtained in (Gurzadyan & Savvidy 1984, 1986) via 
Maupertuis principle. 



2 



V.G.Gurzadyan and A.A.Kocharyan: Relaxation 



is responsible for reaching a coarse-grained state in non- 
stationary systems (violent relaxation) (Lynden-Bell 1967). 

Numerical studies possess their own difficulties, starting 
from the choice of descriptors up to the interpretation of the 
results, e.g., as shown in (Gurzadyan & Kocharyan 1994), 
the computer image of the Lyapunov exponents is non- 
equivalent to their definition. So, the importance of fur- 
ther searches of strict methods to study the chaos in stellar 
systems is doubtless. We revisit the problem of both N- 
body relaxation and its time scale using both geometric 
(Kocharyan 1993) and Van Kampen stochastic differential 
equation approaches (Van Kampen 1976). The results con- 
firm those of the Maupertuis principle for the collective 
relaxation timescale of collisionless spherical systems 3 . 



2. Stochastic instability 

We consider an TV-body system described by the 
Lagrangian 



1 N 

£(x,v) = ~^2 m aWa\ 2 - V(x), 
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where 



f(x) = -G£ 



a<b 



m a m b 
x - Xi,| 



(2) 



(3) 



Hereafter we use units G = 1 and m a = 1, and 
x= (xi,...,xjv), v = (vi, . . . ,vjv), 
are the coordinates and velocities of stars, respectively, 

X a = (x a ,X 2 ,X 3 ), V a = (V^V^V 3 ), 

where a = 1, . . . ,N. 

According to the theory of dynamical systems, the sta- 
tistical properties of the systems can be studied from the 
behaviour of close trajectories (Arnold 1989). It is shown 
(Kocharyan 1993), that the evolution of the distance be- 
tween two nearby trajectories, denoted by £, is described 
by a generalized Jacobi equation, which can be written in 
the form 

i+B(x,v)£ = o, (4) 



where 



B(x,v) 



1 



3N- 1 



Vy^(x) 

2(E-V(x)) 



+ 



3N-1 



AV(x) 

£ l^a(x)| 2 (e^(x)-V c 

a— 1 \a— 1 

2(E - V(x)) A{E-V{x)f 



2 -i 



E is the total energy of the system 



1 N 

E=^l + V(x), 



(5) 



a=l 



For a study of systems with rotational momentum 
(spiral galaxies) with quite different dynamics, see e.g., 
(Gurzadyan & Kocharyan 1988) 



and 



6=1 
b^a 



Xb - X a 
]X6 - xJ J 



(6) 



As in Gurzadyan & Savvidy 1984, 1986, we assume that 
the system is collisionless 4 . We then have AV(x) = and 
one can substitute 



v 4 vi 

v a v b 



to obtain 
where 

W (x) = <B(x,v)> 



2{E-V{x)) W 
£ + lo(x)£ = 0, 



(7) 

(8) 



N 



2N(E-V(x)) 



£|F Q (x)| 2 >0. (9) 



One observes that E — V(x) is the total kinetic energy of 
the system 



1 N 1 

£-v(x) = 2£K| 2 ~-AV>. 



o=l 



Thus, 



JV 



w(x) 



N 2/ y 2\ 

1 ' o=l 



(10) 



(11) 



For spherically symmetric systems, we can replace |-Fa(aO| 2 
with |.F(a:)| 2 (see Chandrasekhar 1942, Cohen 1975), where 



N 



(12) 
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We then replace u> with a stochastic process (cf. 
Casetti et al. 2000, Chandrasekhar 1943). Let X lt ...,X N 
be a sequence of N independent and identically distributed 
(i.i.d.) random variables each having finite values of expec- 
tation \x and variance a 2 > 0. The central limit theorem 
states that as N increases, the sum of N random variables 
given by 

S N =X 1 + ---+X N (13) 

approaches the normal distribution n(/x, a). Then, at large 
N, Sn will behave like a Gaussian process. This can be 
written as 

^#~»(0,1). (14) 

(TV A 

In our case, we have X = |F(x)| 2 , 

M = (X) = (|F(x)| 2 ), a 2 = Var(A) = <|F(x)| 4 ) - M 2 . 

(15) 

Thus, 

w~ w + win(0, 1), (16) 

where 

(17) 



u> 



N(v 2 ) 



2\' 



N 3 / 2 (v 2 )' 



4 By collisionless systems, we understand, as usual, systems in 
which the direct impact of two stars have no role in their dynam- 
ics, and not the neglect of gravitational encounters (scattering) 
of two stars. In this sense, collisionless are even the dense cores 
of star clusters and galaxies. 
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We investigate Eq.(8) by means of a technique 
developed by Van Kampen (Van Kampen 1976) (cf. 
Casetti et al. 2000). One can derive the second moments 
of £ by simply rewriting Eq.(8) in the following form 

H t <* 2 W> \ ( o o 2 \ / (e 2 (t)) \ 

- (l 2 (t)) = 2rul -2^o (l 2 (t)) , (18) 
where (Casetti ct al. 2000) 



1 

T = — 



(19) 



2 2y/u) (u) + UJi) + TTLOl 

The system given by Eq.(18) has a positive Lyapunov ex- 
ponent x defined by 

4w \ 



(20) 



where 



q = ^2fa>2 + v /(2fo;?) 2 + (4^ /3) 3 ) ' . (21) 

We now estimate \x and a to calculate ui n , u)\, and then 
f and x- Since one has (Cohen 1975) 



\F(> 



n(v 2 ), 



(22) 



where n is the mean concentration of stars in the system, 
then 

a 2 = (|F(x)| 4 ) - M 2 ~n 2 N 2 (v 2 ) 2 . (23) 

Therefore, 



1 



UJ-l 



(24) 



And finally we derive the relaxation time 

Tcr-X" 1 -T^N 1 / 3 , (25) 

confirming the result of Eq.(l) derived in 
(Gurzadyan & Savvidy 1984, 1986). 
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3. Conclusion 

The stochastic equation approach used above comple- 
ments the probabilistic approach of Chandrasekhar and von 
Neumann (1943). Both are supported by the fact of decay 
of the time correlation function due to the exponential in- 
stability and the Holtsmark distribution of the fluctuating 
force. 5 

Thus, the stochastic equation method confirms the 
purely geometric derivation of the collective relaxation time 
given by Eq.(l). Although that formula is also supported 
also by alternative numerical analyses (Antoni ct al. 1999), 
the present derivation avoids any approximations in numer- 
ical simulations. 

Chaotic effects could be useful in constraining observa- 
tionally modified gravity theories in the Solar system and 
galaxies (see Capozziello & Tsujikawa 2008). 

5 At distances of the order of the radius of gravitational 
influence r h — l/(v 2 } the Holtsmark distribution diverges 
and should be cut off, as done in (Chandrasekhar 1943, 
Gurzadyan & Savvidy 1984, 1986). For real stellar systems how- 
ever, this cutoff is insensitive to the precise value of rh, since for 
them rh is far smaller than the mean interstellar distance d and 
hence the Holtsmark law vanishes far earlier. 



